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We investigate the relaxation and decay of a particle during inflation by implementing the dynam- 
ical renormalization group. This investigation allows us to give a meaningful definition for the decay 
rate in an expanding universe. As a prelude to a more general scenario, the method is applied here 
to study the decay of a particle in de Sitter inflation via a trilinear coupling to massless conformally 
- - - coupled particles, both for wavelengths much larger and much smaller than the Hubble radius. For 

\l ' superhorizon modes we find that the decay is of the form r;^^ with rj being conformal time and we 

>r^ , give an explicit expression for Fi to leading order in the coupling which has a noteworthy interpre- 

tation in terms of the Hawking temperature of de Sitter space-time. We show that if the mass M 
of the decaying field is <^ H then the decay rate during infiation is enhanced over the Minkowski 
space-time result by a factor 2H/nM . For wavelengths much smaller than the Hubble radius we 
find that the decay law is e'^w'^f''' with C{ri) the scale factor and a determined by the strength of 
the trilinear coupling. In all cases we find a substantial enhancement in the decay law as compared 
to Minkowski space-time. These results suggest potential implications for the spectrum of scalar 
density fluctuations as well as non-gaussianities. 
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I. INTRODUCTION 
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f^ I Inflation was originally proposed to solve several outstanding problems of the standard Big Bang model Q, y, y, 

^^ ' ^ 1^ thus becoming an important paradigm in cosmology. At the same time that inflation solves these problems it 

^£ I also provides a natural mechanism for the generation of scalar density fluctuations that seed large scale structure, 

r^ , thus explaining the origin of the temperature anisotropy in the cosmic microwave background (CMB), as well as 

i' tensor perturbations (primordial gravitational waves). Recently the Wilkinson Microwave Anisotropy Probe (WMAP) 

J^^ ^ collaboration has provided a full-sky map of the temperature fluctuations of the cosmic microwave background (CMB) 

5^ with unprecedented accuracy and an exhaustive analysis of the data confirming the basic and robust predictions of 

^ ; inflationlillQ. 

Cd During inflation quantum vacuum fluctuations are generated with physical wavelengths that grow faster than the 

►^ Hubble radius, when the wavelength of these perturbations crosses the horizon these perturbations freeze out and 

• I— I , decouple[3,0,l3 ■ Wavelengths that are of cosmological relevance today re-enter the horizon during matter domination 

^^ ' when the scalar (curvature) perturbations induce temperature anisotropics that are imprinted on the CMB at the last 

scattering surface 9, lOj. Generic inflationary models predict that these are mainly gaussian adiabatic perturbations 

with a spectrum that is almost scale invariant. 

Inflationary dynamics is typically studied by treating the inflaton as a homogeneous classical scalar field 0, S lj| 
whose evolution is determined by a classical equation of motion, while quantum fluctuations of the scalar field around 
the classical value are treated in the gaussian approximation and provide the seeds for the scalar density perturbations 
of the metric. The quantum field theory interpretation of the classical homogeneous field configuration that drives 
inflation is that it is the expectation value of a quantum field operator in a translational invariant quantum state. 
There are important aspects of the dynamics that require a full quantum treatment for their consistent description, 
for example particle production and in particular particle decay. A systematic treatment of the quantum dynamics 
of the inflaton that includes particle production within a non-perturbative framework is given in ref. pdl ll^ . 

While the dynamics of particle production during inflation has received much attention, the full quantum treatment 
of particle decay during the inflationary (or more generally during a rapidly expanding) stage has not been the focus 
of similar attention. 
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While most reheating mechanisms rely on the coupling of the inflaton to other fields into which it can decay leading 
to a radiation dominated stage, the consequences of such coupling between the inflaton and other fields during the 
inflationary stage are typically neglected. In this article we focus on the study of the decay of a particle that could 
be the inflaton as a result of its coupling to other fields. 

To the best of our knowledge a preliminary study of the decay of the inflaton during a de Sitter stage has been 
previously addressed within a particular case in ref.|l.1|. 

There could be several potentially important consequences of particle decay during inflation: if the inflaton couples 
to other particles, then its quantum fluctuations which seed density perturbations also couple to these other fields. 
Therefore the decay of the quantum fluctuations of the inflaton may result in a modification of the power spectrum 
of density perturbations. Furthermore the coupling of the quantum fluctuations of the inflaton and consequently 
of density perturbations to other fields may possibly induce non-gaussian correlations. For a recent review on non- 
gaussian correlations generated during inflation see ref . [lj| . 

While a thorough assessment of these potentially relevant phenomena requires a detailed treatment of the coupling 
of gauge invariant perturbations to other fields into which these fluctuations can decay, in a spatially flat gauge there 
is a direct relation between the evolution equations for density perturbations and those of the quantum fluctuations 
of the infiaton!l(J|. 

Therefore if density perturbations also couple to these other fields as a consequence of the coupling of the inflaton 
field to these fields and if density perturbations decay into these fields, such decay implies that the amplitude of 
density perturbations will diminish with a consequent modification of their power spectrum. 

Clearly a first step in the program to assess these potential observables is to understand the decay of fluctuations 
during inflation which is the focus of our study in this article. 

These possibilities with distinct potential phenomenological consequences for CMB anisotropics motivate us to 
study in detail particle decay during an inflationary stage, which we take to be described by a de Sitter space-time. 

The decay of the inflaton during a post-inflationary stage has been considered recently[i3 as a possible source of 
metric (and therefore temperature) perturbations arising from an inhomogeneity of the inflaton coupling. However 
most of these treatments rely on the concept of the decay rate of a particle in Minkowski space-time seemingly 
uncritical of its validity in the (rapidly) expanding universe. 

Goals of this article: In this article we study the decay of a particle into other particles during inflation. The 
decaying particle could be the inflaton but our study will be more generally valid. The main focus of our study is to 
provide an understanding of the concept of decay of a particle in a rapidly expanding cosmology, and to introduce 
and implement a method that allows a systematic and unambiguous study of the relaxational dynamics of quantum 
fields and in particular allows to extract the decay law resulting from interactions. In Minkowski space-time there 
are two alternative but equivalent manners to define the decay rate of a particle: I) the total decay rate is the 
inclusive transition probability per unit time from an initial 'in' state to final 'out' states, II) the total decay rate is 
the imaginary part of the space-time Fourier transform of the self-energy of the particle evaluated on the particle's 
mass shell and divided by its mass-shell energy. Both definitions are equivalent by dint of the optical theorem, or 
alternatively, unitarity. The calculation of a total decay rate from definition I) involves calculating the transition 
amplitude from some initial time ti — > — oo to a final time tf — > -|-oo and multiplying by its complex conjugate. In 
Minkowski space-time the transition amplitude from an asymptotic state in the past to an asymptotic state in the 
future is proportional to an energy conserving delta function. In squaring the amplitude, the square of this delta 
function is interpreted as the total time elapsed in the reaction (T) multiplying an energy conserving delta function. 
Dividing by the total time of the reaction [T) one extracts the decay rate. The calculation of the decay rate from 
the total width via definition II) requires that the self-energy be a function of the time difference and invokes energy- 
momentum conservation at each interaction vertex. The space-time Fourier transform of the self-energy features 
branch cut singularities in the complex frequency plane and the imaginary part across these cuts at the position of 
the particle mass shell gives the decay width or decay rate. The important point in this discussion is that in both 
cases the concept of a decay rate relies heavily on energy (and momentum) conservation. Herein lies the conceptual 
difficulty of extrapolating the concept of a decay rate (an inclusive transition probability per unit time) to the case of 
a rapidly expanding cosmology where there is no global timelike Killing vector associated with conservation of energy 
even when there may be space-like Killing vectors associated with spatial translational symmetries and momentum 
conservation. Such is the case for spatially flat Friedmann- Robertson- Walker cosmologies. The manifest lack of energy 
conservation in an expanding cosmology makes possible processes that would be forbidden in a static space-time by 
energy conservation [l^. In addition, contrary to Minkowski spacetime, cosmological modes in general do not decay 
exponentially with time, therefore the definition of the decay rate requires the kind of analysis we provide here. 

The method: 

Particle decay in de Sitter space-time has been previously studied in reference |13| for some very special cases that 
allowed a solution of the equation of motion. In this article we introduce a method that allows to study the relaxation 
of quantum fields and particle decay in great generality. The main strategy is to study the effective equations of 



motion of the expectation values of fields as an initial value problem in linear response including the self-energy 
corrections. The solution of the equations of motion lead to an unambiguous identification of the decay law from the 
relaxation of the amplitude of the field as a consequence of the self-energy corrections (interactions) . When self-energy 
corrections are included the equations of motion become non-local (non-Markovian) and cannot be solved in closed 
form in general. 

When a perturbative solution of the equations of motion is attempted there emerge secular terms, namely terms 
that grow in time and invalidate the perturbative expansion. These secular terms indicate precisely the relaxation (or 
production) time scales. We implement the dynamical renormalization group introduced in 17] to provide a systematic 
resummation of these secular terms leading to the correct description of relaxation and decay. Such program has been 
successfully applied to a wide variety of non-equilibrium situations in Minkowski space time (see jl7| and references 
therein). 

In this article we generalize this approach to study particle decay during inflation. As a prelude to studying more 
general situations, we begin this program by implementing this method to study the decay of a massive and minimally 
coupled particle into conformally coupled massless scalars via a trilinear interaction vertex. After extracting the decay 
law to lowest order in the loop expansion for the self-energy, we study the limit of Minkowski space-time and show 
that the results obtained reproduce those familiar in Minkowski space time. 

Brief summary of results: We introduce the dynamical renormalization group[l3 method to study the relaxation 
of the expectation value of quantum fields as an initial value problem in the general case. 

After introducing the method and discussing its systematic implementation in the general case we illustrate its 
application and study the decay of a massive particle (it could be the inflaton) coupled to conformally coupled 
massless particles via a trilinear vertex in de Sitter space time. This simpler setting allows to present the main 
aspects of the program as well as reveal the important features associated with the expansion in a clear manner. The 
relaxation and decay law is studied to lowest order in the coupling both for wavelengths that are inside and outside 
the Hubble radius during inflation. The decay constant for superhorizon modes have an interesting interpretation in 
terms of the Hawking temperature of de Sitter space-time. In all cases we find that the decay is enhanced during 
inflation as compared to the Minkowsky space-time result. The decay law for modes deep within the horizon feature 
a wavevector dependence that leads to a larger suppression of the amplitude for longer wavelengths. 

The article is organized as follows: in section II we introduce the model and the non-equilibrium Green's and 
correlation functions in arbitrary vacua, which are necessary ingredients for obtaining the effective equations of 
motion. In section III we obtain the equations of motion including self-energy corrections in the loop expansion and 
introduce the dynamical renormalization group method to extract the relaxation and decay law. In section IV we 
study specific cases up to leading order in the interaction and compare to the results in Minkowski space-time. In 
section V we summarize our results and conclusions and discuss potential implications of our results on the power 
spectrum of density fluctuations and non-gaussianity. 

II. THE MODEL 

We consider a spatially flat Friedmann- Robertson- Walker (FRW) cosmological space time with scale factor a{t), in 
comoving coordinates the action is given by 



with 



7^ = 6 



(- + -) 2 



being the Ricci scalar and ^ an arbitrary coupling to the Ricci scalar: i^ = corresponds to minimal coupling and 
^ = 1/6 corresponds to conformal coupling. The linear term in (/i is a counterterm that will be used to cancel the 
tadpole diagram in the equations of motion. 

It is convenient to pass to conformal time rj with drj = dt/a{t) and introduce a conformal rescaling of the fields 

a{t) <j>{x,t) = x{^,v) ; a{t) Lp{x,t) = 5{x,ri) . (3) 

The action becomes (after discarding surface terms that will not change the equations of motion) 
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X,S = d'xdiil- x' -i'^xf-MUv)x^ + S'--iVSf-Ml{v)S^ -gCiv) X S' -C^iv) Jiv) X> , (4) 



with primes denoting derivatives with respect to conformal time 77 and 
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and C{ri) = ait{r])) is the scale factor as a function of conformal time. For inflationary cosmology the scale factor 
describes a de Sitter space-time, namely 
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ait) = e 

with H the Hubble constant and conformal time rj is given by 
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where vo corresponds to the initial time t — 0. We choose 
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During inflation the effective time dependent masses of the fields are given by 
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The Heisenberg equations of motion for the spatial Fourier modes of wavevector k of the fields in the non-interacting 
ig = 0) theory are given by 



where 
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The Heisenberg free field operators can be expanded in terms of the linearly independent solutions of the mode 
equation 



S'l{k-r,)+[e-\(v^--J\^S.ik-ri) = 0. 



Two linearly independent solutions are given by 
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5, (A:; 77) = -i-^-^-^^H^^Hkr^) 
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where Hi ' (z) are Hankel functions. For wavevectors deep inside the Hubble radius —krj >> 1 these functions have 
the asymptotic behavior 



g,(fc;77)''''= °°^e-*'='' 



/,(fc;?7)''''^°°ie*'='' 
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and are normalized so that their Wronskian is given by 

Wig^^ik; v),fuik; v)] = glik; v) luik] v) - 9Ak; v) flik] v) = -i 
The most general solution of the mode equations p2|l is given by 

Suik;v) = Ci(fc;?7o) Uik;v) + C2(fc;??o) 9vik]v) , 



(15) 

(16) 

(17) 



where the coefBcients Ci_2 are determined by an initial condition on the mode functions Su{k\rj) at conformal time 
7?o, namely 

Ci(fc;?7o) = -■i[9,^{k;vo) Sl{k;r]o)- gl{k;r]o) S^{k-r]o)] (18) 

C2{k;r]o) = ~i[fl{k;r]o) S^{k;rio) - f^{k;T]f)) Sl{k;r]o)] . (19) 

The spatial Fourier transform of the free Heisenberg field operators Xfc('7): %('?) are therefore written as 

Xkiv) = a^S,^{k;v)+al^^S:jk;rj) (20) 

hiv) = Pk s.s (fc; '?) + ^U St, (fc; ^) , (21) 

where the Heisenberg operators aj, at and (3^, (3- obey the usual canonical commutation relations. 

Canonical commutation relations both for the Heisenberg fields (and their canonical momenta given by their deriva- 
tives with respect to conformal time) and the creation and annihilation operators entail that 

\C2{k-m)\'-\Ci{k;7jo)\^ = l. (22) 

The vacuum state |0) is annihilated by a^r, Pj:. However a different choice of the coefficients Ci^2 determine a different 
choice of the vacuum state, the Bunch-Davies vacuum corresponds to choosing C2{k;r]o) — 1, Ci(fc;77o) — |18l |. 
Heretofore we generically refer to the different choices of the coefficients Ci_2 constrained by Eq. H22|) . as S- vacua. 
An illuminating representation of these coefficients can be gleaned by computing the expectation value of the number 
operator in the Bunch-Davies vacuum. Consider the expansion of a scalar field either in terms of the Bunch-Davies 
basis gk(j]) ot in terms of the generalized basis S{k]ri), namely 

XkM = H gu^ (fc; ry) + a^- g*^ (fc; ?]) = aj: S^^ {k; tj) + a^- S*^ (fc; ry) , (23) 

the creation and annihilation operators are related by a Bogoliubov transformation 

ai^C2ai-Cia_j: , aj: = C; aj: - C^ aK . (24) 

The Bunch-Davies vacuum |0)b_d is annihilated by a^, hence we find the expectation value 

BD{0\aiaj:\0)BD = l^ip =. N^ . (25) 

Where Nk is interpreted as the number of S-vacuum particles in the Bunch-Davies vacuum. In combination with the 
constraint (|22|) the above result suggests the following illuminating representation for the coefficients Ci^2 



C2{k)^V'^ + Nk ; Ci(fc) = 0Vfee'^'= (26) 

where Nk and 9k are real functions. 

A. Non-equilibrium Green's and correlation functions 

The main ingredients in the program to obtain the decay of the amplitude of the scalar field are the non-equilibrium 
Green's and correlation functions. Consider a generic free scalar field $ quantized with the expansion 

$,,(,7) = a^ S,{k- ,7) + at__ s:{k- rj) . (27) 

The non-equilibrium Green's and correlation functions are given by 

Gt+irhv') = (0|r(<i>,,(r7)<l>_,,(77'))|0)=5.(A;;r7)5:(A;; 770 6(77- 77') + 5:(A;;r7)5.(fc; 77') 6(77'-,;) (28) 

Gnv,v) = (0|f (<i>,,(7;)$_,-(77'))|0)-5.(A;;7;)5:(A;;7;')e(77'-77) + 5:(A;;7;)5.(fc;77')e(77-7;') (29) 

G^+(77,77') - {0\^j:irj)^_^irj')\0) ^ S,{k;rj) St{k;rj') (30) 

Gt-iv,v') - (0|<i>_£(r7')%W|0) - [G,T+(77,77')]* - 5:(fc;77) S^k^r,') , (31) 



where T, T stand for the time and anti-time ordering symbols respectively. These Green's and correlation functions 
are not independent since they fulfill the following identity 

G++(7y,77') + G,T-(r,,77')=G,T+(77,77') + G+-(77,7y'), (32) 

which can be trivially verified. 

For generalized S-vacua, we find for example 



Gt-iv,r]') = [l + Nk]gtik;v)9Ak;vl+Nk9Ak;v)g:ik;i) + ^/[lTNk\Nk [er^'^gAk^v) 9Ak;v') + c.c] . (33) 

For physical wavelengths that are much smaller than the Hubble radius, namely for fcry, krj' 3> 1 the two point 
correlation function above has the following behavior 

G+-(^,ry') ''''t-^^ i±^ e-'=(.-V) + ^ ,^Hv-n') + 1 ^[l + N,]N, cos [^^ + r/) + 0^)] . (34) 

The first two terms are similar to the Wightman function of a scalar field in a bath in equilibrium, whereas the last 
terms with the sum of the conformal times are a distinct feature of the mixing between particle and antiparticle states 
of the Bunch-Davies vacuum. We can impose asymptotically, for wavelengths deep within the Hubble radius that the 
physics be locally that of flat Minkowski space-time with a timelike Killing vector which implies that the two point 
function be translational invariant in time. This condition requires that the occupation numbers fulfill Nk — > as 
k —^ oo. In particular if we further demand that the number of Bunch-Davies particles in a generalized S- vacuum be 
finite, then it must be that Nk < O [l/k^+'') , e > as fc ^ oo. 

III. EQUATIONS OF MOTION AND THE DYNAMICAL RENORMALIZATION GROUP 

As mentioned in the introduction, in Minkowski space-time the decay rate of a particle can be obtained either from 
the imaginary part of the space-time Fourier transform of the retarded self-energy or alternatively from the transition 
probability per unit time. These methods are equivalent by dint of the optical theorem (unitarity) and both rely on 
energy conservation. In a rapidly expanding cosmology energy is not conserved (albeit there is covariant conservation) 
and the self-energy is not a function of the time difference but is a rapidly varying non-local function of the time 
arguments. In this situation a decay 'rate' cannot be identified unambiguously. 

In ref.|l7j we have introduced a method to study relaxation directly in real time. The method relies in obtaining 
the equation of motion of expectation values of the fields in linear response. The decay law (more general than a 
decay rate) is extracted from the time evolution of the am,plitude of the expectation value. The equation of motion 
for the expectation value includes the non-local self-energy contributions in a consistent loop expansion. Although 
the equation is linear, it is generally non-local and when the self-energy is not time translational invariant (no energy 
conservation) it becomes an integro-differential equation which in general cannot be solved in closed form. The method 
introduced in ref.[l7| relies on a perturbative expansion of the solution in terms of the coupling constant. However, 
such an expansion features secular terms, namely terms that grow in time and invalidate the perturbative expansion. 
The dynamical renormalization group ^l^ provides a consistent resummation of the series that leads to a uniform 
asymptotic expansion. The DRG improved solution directly allows to extract the decay law. This method has been 
applied and its applicability and reliability has been tested in a variety of equilibrium and non-equilibrium situations. 
The method is rather general and allows to resum secular terms for any set of differential or integro-differential 
evolution equations. 

While the method has not yet been applied to the case of an expanding cosmology, we will confirm its reliability 
by analyzing the results in the limit when the expansion rate vanishes, namely, Minkowski space-time. 

The method to study decay and relaxation begins by obtaining the equations of motion for the expectation value 
of the decaying field, which is obtained from the non-equilibrium generating functional[llL ll^ 

Z = fvx^ V5^ g.{A[x^5+]-A[x-,^-]} (35) 



Where the superscripts ± refer to the forward (+) and backward (— ) branches along a closed contour in time 
corresponding to the unitary time evolution operator and its inverse, both of which are required to obtain expectation 
values in an initial density matrix (which could describe a pure state) . For details on this method the reader is referred 
toref.[ll[lp. 

Our goal is to obtain the equation of motion for the expectation value of the field x- For tliis purpose, we implement 
the tadpole method by performing the following shift in the spatial Fourier transform of the field y IIILII^ 

xi(r7)=X,~(77)-Ka±(,7) ; {x^{ri)) ^ Xj:{^) ; (aJ(r;))=0. (36) 
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FIG. 1: Self energy loop of conformally coupled massless particles (dashed lines). 



in the above expressions (■ ■ ■ ) stand for expectation values in the initial state which can be prepared by switching 
on an external source term to displace the field and switching the source off to let the field evolve. This is the usual 
method to prepare an initial value problem in linear response. For more details on this method we refer the reader 
to|llul7ij. 

In terms of the shifted field, the action for the spatial Fourier transformed fields becomes 



± x± 



A\X,a^,5- 
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where we have neglected terms that cancel in the difference A[X,a~^,S~^] — A[X,a^ ,d^], and V is the comoving 
spatial volume. The equation of motion for the expectation value Xz{r]) is obtained by implementing the condition 
(cr (77)) = order by order in perturbation theory. The equation of motion is the same for the + and — conditions 
as a consequence of the identity (|32|l . and up to one loop order it is given by 



X'iiv) + [k^+Mliv)] X~^{r,) + 2.g2 ^(77) / dr^' C{r,') ICk{v,v') X^{f^') = , 



(37) 



where we have used the counterterm J{ri) to cancel the tadpole (S^S^g). At one loop order in the fields S {0{g'^)) the 
non-local kernel lCk{rj, rj') is given by the following expression 



!^k{v,v') = -« 



(27r)3 



G^:.*('?'^') ^I^^-.^^^V') Gjziv,v') G^-^.M) 



(38) 



where the correlation functions Gg }~^{r], rj') are those for the 5 field given by the expressions Eas. H28|l - 131|l with v — us 
given in Eq. I|ll|) . This non-local kernel corresponds to the one-loop retarded self energy in real time as shown in fig. 

Ill 

Although the equation of motion (|37|) is linear, it is non-local and a general solution is unavailable. However for 
weak coupling g^ a perturbative solution can be found by writing 



Xiiv) = Y.i9'T x^^M 



n=0 



leading to the hierarchy of coupled equations 






U<-\) 



^'--M-\ 



V 



x,rM = , 

Xr,rM=nn{k-lj) ■ n-1,2. 



7^„(fc;r7) = -2 C(7?) / d^' C{ij') ICk{v,v') ^„_i,^(^') 



(39) 

(40) 
(41) 
(42) 
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In terms of two linearly independent solutions of the unperturbed equation g^, {k,r]);fi, {k,ri) with Wronskian [see 
Eq. (|16l) ] equal to —i, we write 

X^j:{f^) = Aj: g,^ (fc; ry) + i?,, /.^ (fc; ry) , (43) 

where the coefficients A-r^ B^ are determined by initial conditions. The hierarchy of coupled equations H4U|) can be 
solved iteratively by introducing the (retarded) Green's function of the second order differential operator on the left 
hand side of these equations. This Green's function is given by 

g{k; T],r]')^i [g^^ (fc; 77) ^^ (fc; ry') - f^,^ (fc; ry) g^,^ (fc; t]')] 6(77 - 77') . (44) 

Hence the solution of the hierarchy of equations for n > 1 is given by 

K,kiv) = I dv' g{k; V, v') 7^„(fc; v') ■ (45) 

The general form of the perturbativc solution Eq. (|45() combined with the linearity of the equation of motion indicate 
that the full solution is of the form 

Xj:(77) = Aj 5,Jfc; 77) [1 + 52 ^i(fc; Ty) + / i^2(fc; ^) + 0(/)] + Sfc /.,(fc; ^) [1 + 5' ^i(fc; »?) + 5^ ^2(fc; ry) + 0(/)] 

(46) 
where the functions Fn{k,ri); Hn{k]rf) are found iteratively from the procedure described above. If the functions 
Fn{k] 77); Hn{k; rj) remain bounded in the limit 77 — + then the perturbativc expansion provides a convergent uniform 
series. However, in general these functions feature secular terms, namely contributions that diverge in the limit 77 — > 
and invalidate the perturbativc expansion. In order to provide a uniform solution valid at all (conformal) times, these 
secular terms must be resummed via the dynamical renormalization group|l7|. 

A. Dynamical renormalization group (DRG) 

The general solution given by Eq. (|46(l highlights that the perturbativc corrections can be interpreted as a renor- 
malization of the complex amplitudes Aj:; i?^, the dynamical renormalization group provides a systematic manner to 
resum the secular divergences in the perturbativc expansion in terms of renormalization of the amplitudes. 

This resummation program begins by extracting the secular terms of the functions Fn{k;ri), Hn(k]rj) from the 
terms that remain finite and bounded in conformal time. Let us write 

Fn{k;r])^Fn,s{k;v)+Fnj{k;r]) ; -ff„(fc; 77) = 7?„,,(fc; 77) + if„j(fc; 77) , (47) 

where Fn^s{k]v)i Hn^s{k]^) are secular, namely diverge in the limit 77 — > 0, whereas Fnj{k;ri), Hnj{k]r]) remain 
bounded for 77 ^ 0. 

The dynamical renormalization group implements a resummation of the secular terms by introducing an arbitrary 
scale ?7 and a wave function renormalization of the complex amplitudes as follows 17J 

Af: = Aj:i^) Zf{f,) ■ Zf{f,) ^l+g^ zf-{f,) + / ^^^^-(77) + 0{g^) , (48) 

B^^ = i?,^(77 )Z|(7y) ; Z?{f,) = I + g^ z^^if,) + / z^^-(7y) + 0(/) . (49) 

The solution of the equation of motion now becomes 

Xii^) - M^) .9.,(fc; ^) [1 + 9^ {Fi{k; fj) + z^^iv)) + 0(/)] +i?,-(77) /.^ (fc; 7y) [l + g' (H,{k- 77) + zf^{fi)) + 0(/) 

(50) 
The coefficients z U are chosen so that they cancel the secular terms Fns{k; rj), i^„^s(fc; 77) at the point rj ^ fj, namely 

4d^) = -PiAk;fi) ; zfj.{f^) = -H,^s{k;fi) ; etc. (5i) 

Therefore, the solution now becomes 

Xj:ifj) = Aj:if^) 5.,(fc;^) {1 + .9' [F^ik;7^) - FUk;v)] + Oig'')}+Bj:{fj) /.Jfc; 77) {l + g' [H,{k;rj) - H,{k-f,)]+0{g')] 

(52) 



This form of the solution can be written in the more iUuminating manner 



Xj^iv) = ^fc(^) 9u, (fc; ^) 1 1 + 5' /' '^^^i^:'^'^ drj' + 0{g') + non-secularj 

+ %(ry) UJk; r;) (1 + .9^ /' ^^^Ak^ri') ^^, ^ ^^^4^ ^ non-secularl (53) 



k\-i/ J^x\-^ ■// I - ' .V / ^^, 



where the non-secular terms are terms bounded in the limit 77 -^ 0. Ea. H53|l reveals that the solution has been 
improved by choosing the scale 77 arbitrarily close to 77. 

The solution X^{ri) is independent of the arbitrary renornialization scale r), namely 



dXf:iv) 



(54) 



df] 

which leads to the dynamical renornialization group equation\V^. To lowest order the DRG equation is given by 

OA^if,) ,^..(^)^£m(M)+o(/)^0, (55) 



drj df) 

9'Bj:{v) ^^'-f-" +0{g') ^ 0, (56) 






To lowest order in g^ the solution of these dynamical renornialization group equations is given by 

Mv) = Mvo) e9'[^i-('^-^'')-^-=('=^'^°)l (57) 

Bj(ry) = Bj(?7o) e9'[^i-('=^'')-^i.= (^-*)l . (58) 

Since the scale fj is arbitrary, we can now set fj = rj and obtain the renornialization group improved solution 

Xj:{v) = A^^{Tj)g,^{k;^)[l + g^F,j{k,T^)+0{g^)]+B^^{7j)U^{k;Tj)[l + g'H,j{k,7j)+0{g^)] , (59) 

Aj:{ri) = Aj:{7)o) es'I^^-C^^'')-^^-^'^^''")! , (60) 

%(77) = Bj:{t]o) e9'[«i.= ('=^")-f^i'=(''^^'"')] , (61) 

the terms in the brackets in Eq. H59() are truly perturbatively small at all conformal times for weak coupling. That is, 
the dynamical renornialization group produces a perturbative expansion which is uniform in time. The exponential 
factors in the complex amplitudes (|60|I - H61|) will determine the decay of the am plit ude. The reliability and power of 
this method have been tested in many different cases and we refer the reader toT7] for a more thorough discussion. 
We now implement this program in several relevant cases. 

IV. SPECIFIC CASES 

We begin our program by implementing the dynamical renornialization group resummation in a simpler case and 
to lowest order in the coupling, namely 0(17^), with the goal of highlighting the main aspects of the program in a 
simpler setting. 

For this we consider the S field (decay product) to be a massless conformally coupled field in its Bunch-Davies 
vacuum, namely m — 0, £,s — 1/6, C2 = 1 and the decaying field x to be massive and minimally coupled, namely 
^^ = . For the case of massless conformally coupled particles with Bunch-Davies vacuum in the self-energy loop 

S^,{k;v) = ^ e'^'^ (62) 

and the non-local kernel is given by 



d^q 



sm 



{q + \k + q\){v - V') 



1 ., . ,., „ / 1 



l^kiv,v')= TTtU riT-^ ^=-— cos[fc(r7-r7')]P (63) 



(27r)3 2q\k + ^ 87r2 \r] ^ r]' 



where V stands for the principal part. We define the principal part prescription as follows 



V 



1 



77-77 



■n-ri J {rj-T]') +{eri')^ 



1 




1 


1 


2 


77- 


- 1]' + ie rj' rj - 


- rj' ~ ie rj' 







10 



(64) 



This prescription for the principal part regulates the short distance divergence in the operator product expansion 
with a dimensionless infinitesimal quantity e independent of time. This Ty'-dependent point-splitting prescription in 
conformal time correspond to a time-independent point-splitting in comoving time t = —jr \og[—H 77] [see Eq.®]: 



rj — r] 



(77 — 77') + (e r7')2 



= H e 



H t' 



1 



-H(t~t') 



ri_e-H(t-t')l" + , 



(65) 



That is, a time splitting of e/ H between the points t and t' for t ^ t' . This choice of regularization is consistent 
with the short-distance singularities of the operator product expansion in Minkowski space-time, and leads to a 
time-independent mass renormalization. Indeed, time-dependent mass renormalizations are allowed in cosmological 
space-times and they are associated with different regularization prescriptions. (For an analogous discussion using 
the moment cutoff instead of point splitting see sec. Ill of ref. 12]). 

Repeating the calculations that follow but with a (conformal) time independent point-splitting e instead of e 77' as 
conformal time separation in Ea. (|64l) leads to a (conformal) time dependent mass renormalization. While there is no 
unique choice of renormalization prescription, we demand a time independent renormalization of the mass, which is 
achieved by the principal part prescription adopted in Eq. (|64|l . 

Even with the simplification of conformally coupled massless fields in the self-energy loop, the study of the general 
case for arbitrary wavevectors is complicated by the fact that the solutions of the equations of motion at zero order 
are given by Eq. H43() with g^, (k; rf) given by Eq. H13|l and v-^ given by Eq. (|11() for S,x — 0. However, progress can 
be made in the following relevant cases: i) k-q, krj' <C 1 corresponding to wavelengths that are larger than the Hubble 
radius all throughout inflation, which is equivalent to taking fc = 0, ii) fc ^ with |A;77| —f 0, this corresponds to modes 
that cross the horizon during inflation and iii) krj, krj' ^ 1 corresponding to wavelengths that are smaller than the 
Hubble radius all throughout inflation. We study each case separately. 



A. Wavelengths larger than the Hubble radius: k — 



In the case A: = it is convenient to choose the following linearly independent solutions of the unperturbed equations 
of motion 



<?.,(0,r7) = (-77)^+ , /.,(0,77) = (-77)'^- 



with 



1 



/3± = 2 ± '^x ' ^x 



/9 _ M2 
in terms of which we write the solutions for the equation of motion at zeroth order as 



(66) 



(67) 



(68) 



where a and b are constant coefficients. 

The retarded Green's function Ea. H44|l necessary to solve the hierarchy of coupled equations, for A; = is given by 



5(77,77') = -TT^ [{~vf^ {~rff- - (-^7')^+ i-vf-] 6(77-77') = ^ 



2i^ 



2 I/y 



V 



6(77-77'). (69) 



We start by computing 7?.i(0, 77) which from eqs. H42|l . (|64|l and H68|) is given by 



7^1(0,77)=-- 



1 



" dri' a i-v'f+ + b {-if- 



+ {e^-e). 



8 TT^ _ff 2 77 Ljj 77' 1] — 1]' + ie rj' 

Expanding the kernel in a power series of the ratio 77/77' and integrating term by term yields 
1 



7^l(0,77) 



{2tt H 7/)2 



[a (-77)-^+ + 6 (-77)''- ] (log e + 7) + a 



(-77)^+ 7^(1 - /3+) + (-770)^^+ J2 T^ 



fe=i 



77 



k~'f3+ V'7o 



11 



oo 



fc=l 



k-f3^ 



(70) 



where 7 = 0.57721 ... is the Euler-Mascheroni constant and ip{z) stands for the logarithmic derivative of the Gamma 
function. 

Inserting Eas. H69|l and H70|) in Eq. 1451) we find that the first order correction X-^ 5(77) is given by 



X.oiv) = 



1 



2 (2 TT ff)2 i/^ 



i-v) 



/3+ 



-6 {—riY^ { [loge + 7 + '^/'(l — /?-)] log — + non-secular 



[loge + 7 + ?/'(l — P+)] log h non-secular > 

m ) 

V 



(71) 



where the non-secular terms are terms bounded in the limit 77 ^ 0. Therefore, to first order in the coupling we find 
the solution of the equation of motion for the fc = mode to be given by 

^ ' " I 

non-secular 






~ o(o ^m2 [log£ + 7 + ^(/3-)]log^ 



2(2 7rif)2 



[loge + 7 + V'(/3+)] log h non-secular 



'70 



(72) 



where we used that /?+ -I- /3_ = 1 [Ea. (|67|l ]. From this expression we can read off the secular contributions 
i^i,s(0,r;), Hi.s{0,Tj) in Eq. (gZj) for k ^ 0, namely 



HiAQ,v) 



1 



2 (2 TT iJ)2 jy^ 
1 



[loge + 7 + ^(/3_)]ln^ 



0(0 „^^[loge + 7 + ^(/3+)]ln^ 
2 (2 TT Hy v^ ?7o 



(73) 
(74) 



And the dynamical renormalization group resummation analyzed in section (jIII A|) above leads to the following 
resummed solution 



Xq{^) = a(%) {-Vf^ 



[loEe + T + V-O-)] 



[1 -f 0(g2)] -f 6(ryo) (-77)^- 



^[logg + 7 + ^/.O^.)] 



[l + 0(.g2)], (75) 



where the terms in the brackets are perturbatively small [©(g^)] and have a finite limit as 77 -^- 0. The above improved 
solution is uniform for all conformal time, however from this solution it is not clear what is the decay rate since 
the unperturbed solutions {—tj)^^ are multiplied by different functions. In Minkowski space time the decay rate can 
be extracted straightforwardly and unambiguously because it describes in general an exponential relaxation of the 
amplitude that multiplies the oscillatory phases. However, in an expanding cosmology and in particular during a de 
Sitter stage, field modes with wavelengths larger than the Hubble radius do not propagate, they either grow or decay 
as a function of conformal (or comoving) time, thus the concept of the decay rate requires further examination. 

As we explain in the following section, the decay of the amplitude can be separated from a mass renormalization 
in an unambiguous manner. 



Decay rate and mass renormalization: 



The relevant question that we must address is how to recognize a decay of the amplitude from a mass renormalization 
in this expression. We write the mass as the renormalized mass plus mass renormalization counterterms in the usual 
form 



Such renormalization results in a renormalization of v^ and of the exponents (3± , namely 

g2 6Mf 






2 v^.R H^ 
g^ 5Ml 



I "^XR 



9 Mj. 



m 



1 



J--Jj^ , /3±,^=-±7.,, 



(76) 

(77) 
(78) 
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We now insert f3± as given Eq.jTSJ in Ea. (|72|l and to order g^ we find 



^o('7) - « i-vr^ 1 



2 (2 TT iJ)2 l>^ 



[log e + 7 + V(/3- ) + (2 tt)^ (5M?] log — + non-secular I | + 

jyo J/ 



+5 (-r;)'^- ("l + ^ ^^ ^^^^, ^ I [loge + 7 + V^(/3+) + (2 ^)2 W/^] log ^ + non-secular j") 



(79) 



In the expression above and in what follows we have suppressed the subscript R to avoid cluttering of notation, but 
it should be understood that all quantities are renormalized. 

We see that mass renormalization cannot cancel both secular terms, and that the mass renormalization correction 
for the growing and decaying solutions have opposite sign . Therefore, choosing the mass renormalization counterterm 
SM^ to be given by 



SM^ 



1 



f2 7r) 



1 



loge + -/+-m(3+)+iP{p.)] 



(80) 



cancels the logarithmic short distance divergence and leaves a finite contribution that multiplies both solutions equally. 
With this choice of mass counterterm, the solution of the equation of motion to first order becomes 



M^) = « i-vf' 



n 

1 + Fi log h non-secular 

Vo 



6 (-77) 



/3- 



77 

1 -I- Fi log \- non-secular 



(81) 



where 



g tauTTz/^ 
16 n v^ H^ 



(82) 



and we used the relation 



i^{l+A-^[l-^ 



IT tan[7r v] 



We can now apply the DRG which exponentiates the secular terms and gives as the DRG-improved solution after 
mass renormalization 



Mv) 



ri 



a{fio) hvf^ [1 + 0{g')] + 6(770) (-77)^^- [1 + 0(/)] 



(83) 



We now clearly see that a decay rate Fi can be unambiguously identified from the contribution that multiplies both 
solutions, whereas the mass renormalization enters with different signs for each solution. 

Since the amplitudes 0(779), 6(770) obey the dynamical renormalization group equations (|55I56|I . the DRG improved 
solution (|83|l is independent of the scale 770, a change in this scale is compensated by a change in the amplitudes. 

The decay rate in de Sitter space-time can be read off from the above expression since by setting for simplicity 
fjo = m 



g-ri Ht ^ g-rds t^ 



(84) 



and Tds ^ HTi. 

Ea. H82|) can be analytically continued for M > | i? as 



g^ tanh 



Fi = 



7P" 



16nmjM1.9 



(85) 



This result agrees with that obtained by a different method in ref.jl; 
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2. Minkowski space-time limit. 

Minkowski space time is recovered in the limit H ^ 0. In such hmit wc find from eqs. 1)67^ and H85|) 

M 



H~>0 



H 







9 


I&ttMH 



tanh 



kM' 

IT 



H^O 



9 



16ttMH 



Therefore, we find in this hnrit 



tTi 



^ p IGtvM ^ ^ 



(86) 



(87) 



which displays the exponential decay of the amplitude with the correct decay rate Tm = g'^/ilQirM) for long- 



wavelength excitations with mass M in Minkowski space-time. Since —rj 



-Ht 



I H up to an overall normalization of 



the field the dynamical renormalization group improved solution in this limit is given by 



^o(t) 



[A 



JMt 



Be 



-iMt-\ 



(88) 



which is the correct solution for a zero momentum excitation in Minkowski space-time. 

Thus, clearly the dynamical renormalization group provides a consistent resummation program to study relaxation 
in a cosmological setting. 

This comparison highlights that the decay rate in de Sitter space time is related to that of Minkowski space time 
as 



Tds = TiH = T 



M 



M , , 
tan TTi^v 



(89) 



For M ^ H (for the inflaton this corresponds to the slow-roll limit) we find that the decay rate during inflation is 
enhanced as compared to that in Minkowski space-time 



TT M 



>1 



M 



(90) 



The decay rate T^is — HTi with Fi given by the result (|85|l has a noteworthy interpretation in terms of the Hawking 
temperature associated with the horizon in de Sitter space-time Th 



H_ 

2ir- 



181 which can be seen as follows. 



The solution of the free equation of motion for the zero mode given by Eq. H68() in comoving time is given by 



X^{t) = A e*'^+* + B e*' 



UJ± 



± 



9 

4 



(91) 



The homogeneous modes are propagating for M > 37J/2. As compared to the case in Minkowski space-time we 
can identify lo as the complex poles that determine the free field evolution. In terms of these complex poles it is 
straightforward to see that for a propagating mode (M > 3H/2) Yds can be written as follows. 



F 



dS 



g2 l + 2iV(w+) 



167r Re(ci;+) 

with N{uj) the Bose-Einstein distribution function at the Hawking temperature, namely 

1 



N{u;) = 



e^H 



1 



(92) 



(93) 



The expression (|92|l is similar to the expression for the decay rate in Minkowski space time at finite temperature and 
to lowest order in the coupling in terms of the pole frequency for the free field [l9|. 

Thus at least in the case in which there is propagation, namely M > 3iJ/2, the decay rate can be identified as 
being that of Minkowski space time at the Hawking temperature. This noteworthy property of the decay rate for 
superhorizon modes has also been discussed in ref. 13]. 



B. Modes that cross the horizon during inflation k j^ , rj ~* 

For arbitrary k our integral equation (|37|) takes the form 

g \^ 1 n di]' {rj — if) cos k{r] — 7]') 



n^^) 






/ 5 x2 1 PdV 



{rj — rj'Y + (e 77') 



l\2 



x^ii) = 



where we used eas. H63l V H64II . 

To first order in g^ the solution X^ j:{ri) given by Ea. (|45|l becomes 



14 



(94) 



^ , , 1 r rf^' r,, ,, r dv" iv'-v")cosk{r^'-7^") 



(27riJ)2 



77' 



r;" (77' - rj"y + (e 77")^ °''= 



(95) 



where G{k; 77, 77') is given by Eg . (|44|) and for simplicity we consider the solution 

with gi, (k; 77) the mode function with Bmich-Davies initial condition. Ea. H95|) can be written in the following form 



X,r.iv)-Aj"^g^^{k;rj')J{r,,rj'), 
J no '/ 



where 



J{v,v') 



Vv 



drj" (77' -77") cosfc(77'-77") 



87r2ff2 y^, ^" (Y _ ^//)2 _,_ (g ^/) 



Ini 



'Hi,l\kr,)Hi,lHkv") 



(96) 



(97) 



Our goal is to evaluate X^ ^^(77) for 77 ^ .In order to achieve this we need g„^ {k; 77) and the integrand in Eq.JHZI) 
for small arguments: 



9u{k;r])'^^ -—— T 



i T(iy) / 2 I I- ,„ /, 9 9N 



Im 



'Hil\kr,)H(lHkr^') 



■n \k rj 

TT V 



1—7 cot TTV f k rj\ r /OONTl 



1 + 0(772,77'' 



Inserting Eq.lIHHl) in Eq.lEHl and ((nH) yields, 



y 7 ^ ''-o- , 2-2 T{v) . ( 770 



(99) 



where 



77 






1 



1 



(ii 



2 7i y '" ^' Ji (1 + if)i - 1 t''+ 

Carrying out the integrations leads to the following result 



r + (e--e). 



(100) 



rj^ b [ — ) — 7p 



7A(i^+-)+7 + lne 



logH + oM 



(101) 



Notice that this 77 -^ behavior turns out to be fc- independent. This is due to the fact that the term k"^ in Eq.JHU 
becomes negligible compared to the -^ term for 7y — > 0~ after the modes cross the horizon. 

After mass renormalization according to Eq. H8()|l the logarithmic short distance singularity In e is cancelled and we 
find for the mode functions the following result up to order g^, 



^o,.^W + 5'^a('?)''=°"A^ ^ 



r(^) 



TTl'^^a k" 



g tauTTi/j^ 

16 TT ly^ H^ 770 



Iog^ + O(r;o) 



(102) 
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This result features the secular term log 77 which is resummed by implementing the DRG as in sec.IV-A with the 
result, 






Adm) 



Fi 



[1 + (.9^)] 



(103) 



where Fi was defined in Ea. H82() . It is clear from this result that the behaviour for 77 ^ 0~ and k ^0 with \kri\ — > 
is the same as that for the case fc = [see Ea. (|83|l ]. This due to the fact that the physical wavenumbers fcry become 
so small for 77 ^ 0^ that they bear no relevance on the late time dynamics. While this result could be expected on 
physical grounds, it is important to see it emerge from the systematic implementation of the DRG method. 

C. Wavelengths much smaller than the Hubble radius: \krj\ 3> 1: 

Anticipating mass renormalization we write the mass in the equation of motion l|37|l using Eq. H80|l . Furthermore 
for \krj\ 3> 1 corresponding to modes with wavelengths much smaller than the Hubble radius during inflation, the 
hierarchy of equations of motion up to 0{g^) is given by 



X';^{rj) + eX^^f:{T^) - 7^l(fc;77), 



with the inhomogeneity now given by 
7^l(fc;r,): 



where JCkii], r/') is given by Eq.||H2l)- 

The solution of the zeroth order equation is 



'^*^' X^j:irj) - 2 G(j^) C drf C{rj') /C.(r;, V) X^^rj') 



Vo 



0,fe^ 



X^j:irj) = Ak e-^'-" + Bk e*"'" 



(104) 
(105) 

(106) 
(107) 



The counterterm 6Mf is chosen to cancel the short distance divergence proportional to In e/rj'^. After straightforward 
but lengthy algebra we find the following expression for the inhomogeneity in the limit |fc77o| 3> \kri\ ^ 1 



7^l(fc;77) = 



1 



87r2772ij2 



A, e-^'^" 



,77 . vr' 
In -L + i 


+ Bk e''^'" 


W-- 


TT 
-I - 


L m 2j 




l Va 


2j 



(108) 



where the dots stand for terms that are subleading in the limit \kri\ ^ 1. The inhomogencous equation for the first 
order correction can now be solved in terms of the retarded Green's function 



QiVi V) = j: sin[k{r] - ry')] 0{r] - 77') . 
To leading order in the limit \krio\ 3> \kri\ 3> 1 we find 

Ake-^( 2 7 77 \ Bk e'^''^ f ^, , r^/ ^ ^ 2 7 77 

^•'^ i'l ^ i-i 1^ ^ n H rj rjQ ) 32 n H k [ n H r] r]o 



32tt H k 



(109) 



(110) 



where again the dots stand for terms that are subleading in the \krio\ ^ jfcT;! ^ 1 limit and C{ri) = —l/Hr] is the 
scale factor. The term C{ri) — C{r]Q) in the above expression is truly a secular term, since it grows by a factor larger 
than e^^ during inflation. The validity of the perturbative expansion for this term is determined by the requirement 
that \k/HC{i])\ = \kri\ ^ 1, namely that the wavelengths are much smaller than the Hubble radius all throughout 
inflation. 

Thus the solution up to this order is given by 



Xj^irj) = Ak e-''^" \ 1 
+Bk e^'^" 1 1 



^ [C(77)-C(77o)]+7 ^ 



32nHk 



^' [C{rj) - CM] - * ^' 



16tt^H^ ki] 
In^ 1 



32TrHk 



IGiT^H^ kri 



(111) 
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where the dots stand for terms that are of higher order in g^ and subleading in the hniit |fc?7o| ^ \kri\ ^ 1. The 
dynamical renormahzation group resummation H59I61|I leads to the following DRG improved solution 

Xj(r;) = e^;524r['='('')-'='(''«)lJ Afee~'['=''+'^'=('')1 [l + C'(/)] + Sfc e''['=''+'^''(")l [1 + 0(5^)] I , (112) 

where fkifl) is a- logarithmic phase that is not relevant for the decay of the amplitude, and the terms in the brackets 
are truly perturbative in the long time limit for wavelengths much smaller than the Hubble radius. In (|112|l we have 
chosen the renormahzation scale 770 to coincide with rjQ. The DRG improved solution (|112|l reveals the decay of the 
amplitude with the scale factor. The result above has the correct limit in Minkowski space-time as can be seen from 
the following argument. In comoving time, the difference C{r]) — C{r]Q) = e^* therefore in the limit H ^ 

^' [C(r;)-CM]''=":^t, (113) 



327rilfc ' ^" " ' " 327rfc 

which gives the correct exponential relaxation of the amplitude of the field for large momentum in Minkowski space- 
time as shown in the appendix. 

The results for the decay laws reproduce the decay rates in Minkowski space time in the limit H ^ [see Eqs. H87|l 
and H113|) ] thus confirming the reliability of the dynamical renormahzation group approach. 

We can summarize the results obtained above as follows. Consider the solution X„ zi'']) ~ gu{k\ rj) of the unperturbed 
equation with Bunch-Davies initial conditions as given in eq. H13|) . The asymptotic behavior of the power spectrum 
(here we do not include the k^ normalization) of the unperturbed solution for modes deep inside the horizon \kri\ 3> 1 
and those well outside the horizon jfcr;! ^ is given by 

.,2lH-o 2^-^r^(t.)7; 
I^O.SWI - . (fc,)2. (114) 

Particle decay modifies the amplitude of the solution and consequently the power spectrum, which after the DRG 
resummation is given by 



2iy-2 7^2/' N r __ n ZFi 



|x.(,)P""^° ^VtSf^4(*) 



(116) 



where we have normalized the mode functions to Bunch-Davics initial conditions at the beginning of infiation t] = tjq 
in Eq. (|115l) . The solution for wavelengths larger than the Hubble radius is independent of the scale ijo because the 
amplitude Aj:{fiQ) obeys the DRG equation Ea. H55() . This amplitude at a given scale 770 is obtained by matching the 
asymptotic forms of the DRG improved solution at a scale fjo. Clearly this amplitude will depend on the decay law 
of modes deep inside the horizon, which reflects a larger suppression of the amplitude for long wavelength modes. 

These results are general hence they are also valid for the decay of the quantum fluctuations of the inflaton fleld. 
Since the quantum fluctuations of the inflaton field seed scalar density perturbations the result obtained above leads 
us to conjecture that the process of particle decay can lead to modifications of the power spectrum of superhorizon 
density perturbations which is obtained when the fluctuation freezes as 77 ^ 0~|lO||. The new renormahzation scale 
770 will lead to violations of scale invariance much in the same way as in the renormahzation group approach to deep 
inelastic scattering. 

Clearly in order to assess the possibility of corrections to the power spectrum of density perturbations from decay 
of quantum fluctuations, the following issues must be studied further: i) a full gauge invariant treatment of the self- 
energy corrections to the equations of motion for density perturbation, ii) a DRG-improved solution for the whole 
range of momenta. Such program is necessarily beyond the scope of this article but the results above are suggestive 
of potentially important corrections to the power spectrum resulting from the decay of quantum fluctuations. 

V. CONCLUSIONS AND DISCUSSION 

The main goals of this article are a study of particle decay in inflationary cosmology, and to introduce and implement 
a method based on the dynamical renormahzation group that allows to systematically obtain the relaxational dynamics 
of quantum fields and the decay law in particular. 
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One of the main points of this work is that during inflation or more generally during a period of very rapid 
cosmological expansion, the concept of a decay rate is ill suited to describe the relaxational dynamics or particle 
decay. In these cases of relevance in Early Universe cosmology, the lack of a global time-like Killing vector prevents 
the interpretation of a decay rate as an inclusive transition probability between asymptotic in and out states per unit 
time and deems unreliable the Minkowski space-time decay rate to describe particle decay. 

The method that we propose to study the relaxational dynamics and to extract the decay law relies on the full 
quantum equation of motion of the expectation value of the field in an initial state in linear response. The quan- 
tum equations of motion are non-local as a consequence of loop corrections which determine the self-energy. The 
perturbative solution of these non-local equations of motion feature secular terms, namely terms tha t gr ow in time 
and lead to a breakdown of the naive perturbative expansion. The dynamical renormalization group I17| provides a 
systematic resummation of the perturbative expansion that leads to an unambiguous understanding of the decay law. 
The dynamical renormalizati on g roup program has been successfully implemented and tested in a variety of situations 
in Minkowski space-time (see \V^) and this work extends it to the case of expanding cosmology, in particular de Sitter 
space-time. 

After introducing the method within a familiar model of interacting fields in inflationary cosmology we studied the 
relaxational dynamics of a massive field whose quanta decay into massless conformally coupled particles via a trilinear 
coupling. This model allows us to present the method and highlight several important aspects in a simpler setting. 
We have studied the relaxational dynamics and the decay law in the following cases: i) fc = 0, namely supcrhorizon 
modes, ii) fixed k ^ and 77 — > 0~ {\kri\ -^ 0), namely modes that cross the horizon during inflation and iii) modes 
deep within the Hubble radius during inflation \kri\ >> 1. 

Cases i) and ii) are found to be equivalent insofar as their relaxational dynamics. The decay constant in this case 
has a noteworthy interpretation in terms of the Hawking temperature of de Sitter space time and the Minkowski 
space-time limit H ^ reproduces the familiar decay rate of a massive particle into massless ones. In the case of 

modes that are deep inside the Hubble radius throughout inflation, we find that the relaxation is exponential in the 

2 

scale factor, the amplitude decays as [see Eq. H112(l ] e~32^^fre[C('?)-C'(''o)l. We have confirmed that the limit iJ ^ 
reproduces the Minkowski space-time result. 

In all cases studied here we find that the expansion enhances the decay. In the case of superhorizon modes we find 
that for H » M (with M being the mass of the decaying field) the rate constant in de Sitter space time is larger 
than that in Minkowski space time by a factor ~ H/M . 

Our results summarized in Eqs. I|115l) - (|116|l for the decay law of modes deep within the horizon as well as those 
that are superhorizon during inflation lead us to suggest potential observational implications. In an interacting 
theory the quantum fluctuations of the inflaton field will decay as a consequence of the interaction. These quantum 
fluctuations (in a suitable gauge ^. IQj) are the quantum seeds of metric perturbations. If these fluctuations decay 
as a consequence of the coupling between the inflaton and other fields (such a coupling is typically assumed for a 
post-inflationary reheating stage) then the decay of the amplitude both for modes deep within the horizon as well 
as those that cross the horizon during inflation will result in potential corrections to the power spectrum of density 
perturbations. In particular the amplitude of superhorizon modes will depend on the decay law of modes inside the 
horizon, which displays a larger suppression of the amplitude for modes of longer wavelength. 

Furthermore, non-linear interactions are necessarily the source of non-gaussian correlation functions, in the case of 
a cubic vertex the three point function in the Born approximation also reveals the emergence of terms that grow in 
time |20j . This three point correlation function is a measure of non-gaussianity, and is clearly of interest to study 
further if and how the decay of inflaton fluctuations leads to non-gaussian correlations, perhaps by implementing the 
DRG as a resummation of the secular terms. 

The decay of the inflaton field into other fields as well as the decay into its own quanta and the implications for 
the density fluctuations will be explored in a forthcoming article [2]| . 
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APPENDIX A: DECAY RATE IN MINKOWSKI SPACE-TIME: 

In order to compare the results obtained for the decay in inflationary cosmology with those more familiar in 
Minkowski space time we now summarize the Minkowski case. The equations of motion in Minkowski space time are 



given 



bylia 



X^{t) + ul X^{t) + / dt' Kk{t - t') X^(t') = 



<-k' 



Ar 



where the non-local kernel is given by 



f rfin sill 



(q+\k + q\){t-t') 



0^ 2q\k + q] 

This kernel can be written in terms of the spectral density crfe(w) in the form 

/+00 
-CXD 

with the spectral density given by 



akiuj) =5^ 



d\ 



16 t:^ q \k + q\ 



S{lu - q - \k + q\) - S{lu + q + \k + q\) 
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(Al) 



(A2) 



(A3) 



(A4) 



The delta functions represent the kinematic cuts for particles and antiparticles respectively. An analysis of the self- 
energy of the decaying particle reveals that the spectral function is related to the imaginary part of the retarded 
self-energy as follows [ 



1 



afc(aj) = — lmS{k;Lu) 



(A5) 



where Iml](fc;cj) is the imaginary part of the retarded self-energy. A straightforward calculation gives the following 
result 



c^fc(w) ^ -^ sign(w) Q{\uj\-k) 



The decay rate F^ is given by 



ImS(fc;wfc) 

2wfc 

where LUk is the particle pole (mass shell dispersion relation) and leads to the result 

o2 



rfc = 



9 



16 TT UJk 



8(wfc - k) 



The zero momentum limit gives 



IGttM 



(A6) 



(A7) 



(A8) 



(A9) 



which is the result quoted in Eq. H87|) for the limit H ^ 0. 

The limit M = which describes the Minkowski space-time limit of the case of wavelengths much smaller than the 
Hubble radius [see section ITVC 'I and Eq. ()104|l ] yields, 

(AlO) 



327rfc 



since 9(0) — 1/2. This result coincides with the i? — > limit in Eq. H113(l . A detailed study using the dynamical 
renormalization grouD|l7j reveals that if the particle mass shell coincides with the origin of the multiparticle threshold 
there emerges a phase that depends logarithmically on time. Such is the case in the massless limit and the logarithmic 
phase is the Minkowski space-time limit of the correction cpk (rj) in the improved solution (|112(l (for more details see 
ref.[l3). 
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